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We study a Luttinger Liquid in a finite one-dimensional wire with box-like boundary conditions 
by considering the local distribution of the single particle spectral weight. This corresponds to 
the experimental probability of extracting a single electron at a given place and energy, which can 
be interpreted as the square of an electron wave-function. For the non-interacting case, this is 
given by a standing wave at the Fermi wave-vector. In the presence of interactions, however, the 
wave-functions obtain additional structure with a sharp depletion near the edges and modulations 
throughout the wire. In the spinful case, these modulations correspond to the separate spin- and 
charge-like excitations in the system. 
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The problem of a particle in a one-dimensional box is a classic example in almost any quantum mechanics text-book 
since it gives a pedagogical introduction to the concept of energy quantization and provides a complete visualization 
of the corresponding wave-functions. It is however only recently that this problem has gained true experimental 
relevance due to the progress in constructing smaller and more refined structures to confine electrons in dots and 
wires. It is for example now possible to resolve the electron wave functions in a finite piece of carbon nanotube 
by Scanning Tunneling Microscopy (STM) experiments Q, |^ Q - Most experimental realizations of one-dimensional 
electron boxes contain many electrons in a Fermi sea, but it is possible to study a single particle excitation on top of 
such a ground state configuration and classify the possible energy levels and wave-functions. 

However, electron-electron interactions may produce interesting effects in such one-dimensional many-body systems 
and systematically change the shape of the electron wave-functions. In faet single particle excitations are no longer the 
eigenstates of an interacting Hamiltonian, but of course it is still interesting to determine the probability of extracting 
or inserting individual electrons at a given position and energy. This local probability density can be interpreted as 
the square of the electron wave-function. We therefore study the fundamental problem of single particle excitations 
in a many body interacting Fermion system confined to a one-dimensional box using the Luttinger Liquid formalism. 

We find that the classic example of a box provides again a good visualization of the effect of interactions on the wave- 
function and the energy quantization. In particular, in addition to the expected rapid Fermi wave-vector oscillations 
in the wave-function we can recognize long wavelength modulations, which correspond to the underlying boson-like 
excitations in the Luttinger Liquid. For repulsive interactions the wave-functions are sharply depleted at the edges 
with a characteristic power-law. Analytic expressions for the wave-functions of the first few levels are presented. 

The Luttinger Liquid formalism is a well-established tool to describe interacting electrons confined to one dimension 
^ 1^ . In a linearized region around the Fermi points the Fermion field can be expanded in terms of left- and right- 
movers 

^{x, t) « e^'^-^^flXa;, t) + e-'-'^-^i^Lix, t). (1) 

The Fourier modes of the loft- and right-moving Fermion density are then represented by bosonic creation and 
annihilation operators, which effectively aet by "shifting" Fermions m steps up or down the Spectrum. In the presence 
of interactions it is then possible to solve the model by a Bogoliubov transformation which mixes the left- and right- 
moving bosons. This transformation can be described by a single "Luttinger Liquid parameter" g, which gives the 
hyperbolic cosine a and sine (3 of the Bogoliubov transformation 
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Here g < 1 for repulsive backscattering interactions and (7=1 for no interactions. Forward scattering does not affect 
this parameter, but rescales the effective Fermi- velocity. 

Let us now consider spinless Fermions in a one-dimensional box of length L with fixed boundary conditions ^'(0) — 
^(L) = 0. After bosonization the Fermion helds become exponentials of the boson operators in a linearized region 
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where a is a short-distance cutoff parameter. The fixed boundary condition therefore relates the left- and right- moving 
boson fieids 4)R{x,t) = —4>L{—x,t) + (/)n and determines the mode expansion in terms of ordinary boson creation and 
annihilation operators and zero modes |i 
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where [4>o,Q] = i and øo = "2^- The eigenvalues of the zero modes Q 



n — riQ — i)y^/(7 are quantized with the 
number of electrons n, where uq = corresponds to the number of electrons in the ground state. 

The Hamiltonian is given hy H = ^ J2m ''^^In^m + ^Q^: where v is the renormahzed Fermi velocity. We see that 
the last term in H resembles a "charging" energy proportional to the square of the excess number of Fermions, but 
this will not afi'ect our calculations since we always consider single particle excitations with exactly one additional 
Fermion n — uq + 1. In general there may also be an additional capacitative energy with a corresponding single 
particle charging energy Eq. 

The boson excitations become highly degenerate with increasing energy levels which are always quantized ujm = 
771^. However, we are interested in the corresponding Fermion wave-function of a single particle excitation on the 
ground state (aj„j|a;) = (aj„j|\E'^^(a;)|0). The probability density p{ujjn,x) is given by the sum of the corresponding 
degenerate wave-functions (ujmlx) squared 



piuj^,x) = V |(tJ„„A|*^(a;)|0)| 
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This is the local density of states which gives the experimental probability of tunneling an electron into the system 
at energy ujm and position x. This spectral density can be readily evaluated for an equally-spaced Spectrum by the 
Fourier transformation of the Fermion Green's function 
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After defining a "mixed wave" Xm {x) = ae*™ l + /3e i , we find 



{ip^{x,t)'ipl{x,0)) = c(sm^'j 
(V'«(x,i)^l(x,0)) . 



2q/3 



exp 



e i- 



Lfc=l 



C I sm I exp 



.fc=i 



-xUx) 



(8) 
(9) 



Here c — ^-j^ i^ir) ^ non- universal cutoff-dependent renormalization parameter which sets the units in our 

calculations and suppresses the spectral weight in the presence of interactions. These correlation functions can be 
simplified to powerlaws of sine-functions ØiØØ but the form above avoids any singularities in the integral of Eq. Q 
since for a given level uj,n we can truncate the sum in the exponential hy k < m. For the first few levels we are even 
able to evaluate p analytically 
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p{u;s,x) = c(sin^)'"'' [x?(x)e''=--] +W [xi{x)x2(.x)e"^-'=] + [x3ix)e^''-^] 
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FIG. 1: The wave-function squared corresponding to m = 1 and ng — 40 in units of c. The modulations become more 
pronounced with increasing interaction strength (smaller g). Last panel: the Integrated probability density for the first 50 
levels compared to the powerlaw behavior ui'^'^ . 



The probability density shows an oscillation of 2kpx with modulations according to the mixing of left- and right- 
moving components in Xm- As can be seen in Fig.^ the amplitude of the modulations increases with the interaction 
strength (smaller g) and the envelope shows a depletion near the edges with a characteristic powerlaw 
agreement with the notion of a boundary exponent 0, IM Ø- limit of small level spacing L — oo we 

recover the known Luttinger Liquid powerlaw behavior of the Integrated density of states p{uj) oc (last panel of 
Fig-ffll- Fig. |21we see that the modulations of level ujm always have m "nodes" and m + 1 maxima with roughly 
equal spacing and height resembling a standing wave with a small wave-vector uj„i/v, corresponding to the density 
waves from the boson excitations relative to the Fermi-energy. It is also instructive to consider the non-interacting 
limit g = 1, for which we always recover a normalized standing wave of wave-vector 2(fci? + LUm/v) without any 
modulations or depletion, even though a general single Fermion excitation still corresponds to a superposition of 
degenerate many-boson states (and vice versa) . 

It is now interesting to explore this modulation pattern for the spinful case where we expect separate spin- and 
charge-like excitations. In this case the electron field with spin a = ± is expressed in terms of spin and charge boson 
operators with two Luttinger Liquid parameters gs and gc 



''pRaix) = giV2TT(a^rj>R,c-Øact>L.c) ^ii7V2TT(as<pB..s-l3s(t>L,B) ^^^^ 

and the analogous expression for 'ijjL,a- The mode expansions for the spin and charge bosons are the same as in Eq. jSJl, 

and the Hamiltonian is also given by the a simple sum H = X)j/=c s (X^m ■^^'^i''^m„'^m„ + ^^l*" ) • Therefore, the 

spin and charge excitations are decoupled except for the quantization conditions on the zero modes Qs = I, 

— {n — \ — ^^^) that [/, n] must be either both even or both odd ([1,1] in our case). For the non-interacting 
case, spin- and charge-excitations are exactly degenerate, but now the excitations are classified by the product space of 
two evenly spaced boson spectra with different energy spacing Vg < v^. The partition function and the electron Greens 
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FIG. 2: The wave-function squared of the first few levels for g — 0.7 and no — 40 in units of c. 



function factorize. Therefore, the wave-functions are products of spin- and charge-modulations which are similar to 
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the ones shown in Fig.Øand Eq. fTn| . except that the mixed wave is rescaled Xm,v{x) = -^e*™ 

and the overall factor is given 

In Fig. Ø we show the square of the wave-functions at the lowest energies in a spinful interacting electron system 
with Qs = 0.7, gc = 0.6, and Vc = 2.3vs- Due to the different velocities, the degeneracy is lifted and many more levels 
appear as the energy is increased [see also Fig. 4 in Ref. 0]. Each level is classified by a spin and a charge quantum 
number and rric, which is reflected by the corresponding number of nodes and maxima in the wave-functions. For 
example the level iris = 4, rric = 1 shows a superposition of two charge maxima and five weaker spin maxima. In Fig.Ø 
we have chosen values of gs = 0.7 and gc = 0.6, which emphasizes the locations of both spin and charge modulations, 
but in systems where the interactions are mostly spin-independent gs is likely to be much closer to unity and the spin 
modulations are less pronounced. Nonetheless, the level structure due to Vs < Vc as well as the charge modulations 
are likely to be observable when classifying Luttinger Liquid systems in real space with STM experiments. 

The Integrated weight of the individual levels decreases with increasing energy, but the total (averaged) density 
of states increases with the known powerlaw p{uj) oc oj^^''^^''. The superstructures survive even in the continuum 
limit and give rise to the observed slow oscillations with wave-vectors lv/vc and lu/vs near the edge of a semi-infinite 
Luttinger Liquid 0, IM Hol llll | . In the non-interacting limit = 5c = 1 we recover again standing waves with a fixed 
wave-vector without any modulations, and remarkably all the many degenerate spin and charge boson states at each 
energy level exactly sum up to an Integrated normalized spectral weight of unity. 

In conclusion we have shown that the single particle wave-functions of interacting Fermions in a box can be used 
to visualize the true nature of the underlying boson excitations. Given the rapid progress in STM imaging and 
nano-structured materials this may soon be observable when classifying different kinds of potential Luttinger Liquid 
systems. 
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FIG. 3: The single particle density of a spinful interacting electron system as a function of energy and position along a finite 
box for Qs = 0.7, Qc = 0.6, and Vc = 2.3'i;a (slightly broadened, so that the fast oscillations are averaged out). This may 
resemble the outcome of a dl /dV measurement in an STM experiment on a finite interacting wire. Top: fuUy resolved level for 
nis ~ 4, rric — 1, and — 40. The two charge peaks and the five weaker spin peaks are visible. 
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